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ABSTRACT 



The dynamical systems of planet-belt interaction are studied by the 
fixed-point analysis and the bifurcation of solutions on the parameter space 
is discussed. For most cases, our analytical and numerical results show that 
the locations of fixed points are determined by the parameters and these fixed 
points are either structurally stable or unstable. In addition to that, there are 
two special fixed points: the one on the inner edge of the belt is asymptotically 
stable and the one on the outer edge of the belt is unstable. This is consistent 
with the observational picture of Asteroid Belt between the Mars and Jupiter: 
the Mars is moving stablely close to the inner edge but the Jupiter is quite far 
from the outer edge. 
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1. Introduction 

The discovered number of extra-solar planets is increasing dramatically due to 
astronomers' observational effort, therefore the dynamical study in this field is getting 
important. Because the belts of planetesimals often exist among planets within a planetary 
system as we have in the Solar System, it is indeed important to understand the solutions 
of dynamical systems of planet-belt interaction. Jiang & Ip [2001] predicted that the 
interaction with the belt or disc might bring the planetary system of upsilon Andromedae 
to the current orbital configuration. 

Yeh & Jiang [2001] used phase-plane analysis to study the orbital migration problem of 
scattered planets. They completely classify the parameter space and solutions and conclude 
that the eccentricity always increases if the planet, which moves on circular orbit initially, 
is scattered to migrate outward. These analytical results is consistent with the numerical 
simulations in Thommes, Duncan & Levison [1999]. 

In addition to astronomy, general or Newton's dynamical systems are studied in 
many other fields and have very important applications. Clausen et al. [1998] studied 
periodic modes of motion of a few body system of magnetic holes both experimentally and 
numerically. Kaulakys et al. [1999] showed that a systems of many bodies moving with 
friction can experience a transition to chaotic behavior. 

On the other hand, Chan et al. [2001] studied bifurcation for limit cycles of quadratic 
systems interestingly. Similar type of approach should be also good for the bifurcation of 
solutions for dynamical systems of planet-belt interaction. In this paper, we focus on the 
planet-belt interaction and study the bifurcation of such system by phase plane analysis. 

Basicly, we would like to understand the orbital evolution of a planet which moves 
around a central star and interacts with a belt. The belt is a annulus with inner radius r\ 
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and outer radius r 2 , where r x and r 2 are assumed to be constants. We set T\ — 3 and r 2 = 6 
for all numerical results in this paper. 

We assume the distance between the central star and the planet is r, where r is a 
function of time. When r < r 1; the belt would only give the force which pulls the planet 
away from the central star. When r > r 2 , the belt would only give the force which pushes 
the planet towards the central star. When r\ < r < r 2 , in addition to the usual gravitational 
force, there is friction between the planet and the belt. These three cases will be studied in 
Model A (r < ri), Model B(r > r 2 ) and Model 0>{j\ < r < r 2 ) individually. 

We will mention our basic governing equations in Section 2. In Section 3-5, we study 
the locations and stabilities of fixed points by phase plane analysis for Model A, B and C. 
The conclusions will be in Section 6. 



The Model 



In general, the equation of motion of the planet is (Goldstein 1980) 



fu _ mf (l) 
d6 2 U l 2 u 



where u — 1/r, to, I are the mass and the angular momentum of the planet and / is the 
total force acting on the planet. We use the polar coordinate (r, 9) to describe the location 
of the planet. 

The total force / includes the contribution from the central star and the belt. The 
force from the central star is 

/. - -% P) 

where we have set the mass of the central star to be one. If the density profile of belt is 
S(r) = co/r where cq is a constant completely determined by the total mass of the belt. 
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The total mass of the belt is 

M m = r P £(r')r'dr'# = 2nc (r 2 - n) (3) 

J Jri 

In general, the force from the belt for the planet is complicated and involved the 
Elliptic Integral. We use a simpler integral to over-estimate the force and then use a small 
number < (3 < 1 to correct the force approximately. 

When the planet is not within the belt region, the force from the belt is 



\ 

-7r Jr 



GmE(r')r' 



-7T Jri r 2 + (r') 2 — 2rr' cos (j) 
r\ — r)(r + r 2 ) 



dr'dcj) 



r \ 11 

ixGcm ( 
= ± < In 



r 2 — r)(r + r\ 
(ri - r)(r + r 2 ) 
(r 2 - r)(r + r x ) 



(4) 



where + stands for the case when r < r\ (the belt pulls the planet away from the central 
star), — stands for the case when r > r 2 (the belt pushes the planet towards the central 
star) and we define c = /3c . 

When the planet is within the belt region, the gravitational force from the belt is 

^ = -^-J. r2 + (r0 2 _ 2^W ^ d<P + (3 L L r* + (rO 2 - 2rr-cos0 rfr 
7rGcm , f (r + ri) (r + r 2 ) 



In 



(r - n)(r - r 2 ) 



4r 2 — e 2 



(5) 



where e is a small number to avoid the singularity of the integral. To be convenient, we 
define the effective mass of the belt to give force on the planet to be 



M b = pM b0 = 27rc(r 2 -n). 



(6) 



Because there might be some scattering between the planetesimals in the belt and 
the planet, we should include frictional force in this case. This frictional force should be 
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proportional to the surface density of the belt at the location of the planet and the velocity 
of the planet. However, if the planet is doing circular motion, the probability of close 
encounter between the planet and the planetesimals in the belt is very small and can be 
neglect here. Thus, we can assume that the frictional force is proportional to the radial 
velocity of the planet dr/dt only and ignore the d6/dt dependence. 

Therefore, f a is proportional to the surface density of the belt and radial velocity 
dr/dt. Hence, we write down the formula for frictional force as 



dr acdr 

fa = ~ aj:{r) d't = 'virv (7) 

where a is a frictional parameter. Because the d6/dt component of the planetary velocity is 
ignored in the frictional force, the angular momentum I is conserved here, so we have 

mr 2 d9 = Idt. (8) 

Because of this, we can use 6 as our independent variable. We use 6 to label time t afterward 
and one can easily gets t from the above equation. 



Since u—l/r and equation (|j), we have 

dr dr d6 I dr I du 



dt d6 dt mr 2 d6 m d9 

Equation (0) becomes 



(9) 



lacudu . . 



3. Model A 



In this model, we assume that the star-planet distance is less than the inner radius of 
the belt, i.e. r < r±. Thus, / = f s + /& + , where f s and fb+ are defined in Section 2. Since 
r < r\ < r%, from equation (|4]), we have 
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. ixGcm \ 
fb+ = < In 



(ri - r)(r + r 2 ) 
(r 2 - r)(r + r x ) 



We further define that u = 1/r and also 

(r x it -!)(! + r 2 w) 



(r 2 M -!)(! + riu) 



We can then transform equation (|I]) into 



(i 2 w _ 

d6t = ~ U + ~P 



Gm 2 iiGcm 2 



{mP/(u)}. 



(11) 



(12) 



Therefore, the equation of motion for this system can be written as 



du 

~d0 

du 

de 



-u + k 2 -—ln(P I (u)), 
u 



(13) 
(14) 



where k 2 = Gm /{I ) and k 3 = (cnGm ) 



3.1. Fixed Points for Model A 



The fixed points (u, v) of problem (0) satisfy the following equations 



and — u + k 2 ln(Pj(w)) = 0. 

u 



(15) 



Obviously, only those solutions locating in the region that u > 1/ri have physical 
meaning since we consider the case that r < r± is this model. However, the properties of the 
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fixed points in unphysical region would affect the topology of the solution curves in physical 
region, so we have to study fixed points both in physical and unphysical regions. 

If we define 

{ u 2 k 1 

Q A (u) = exp j-— + j-u\ , (16) 

from equation (|l5l), we know that fixed points should satisfy Pi{u) = Qa{u). In the 
following theorem, we prove that there are at least three fixed points in whole phase space 
and at least one locates within the physical region, u > l/r±. 

In our theorems, we denote (l/ri) + to represent that u tends to l/r\ from the right 
hand side and (l/ri)~ to represent that u tends to 1/ri from the left hand side etc. 

Theorem 3.1 

(a) There is at least one < l/r 2 such that P/(m*) = Qa{ u *)- 

(b) There is at least one w** G (l/r 2 , 1/ri) such that Pi{u^) = Qa{u^). 

(c) There is at least one u*„ > 1/ri such that Pj(u^) = Qa( u ***)- 
Proof : 

Let Ra(u) = Pi( u ) — Qa(u), where Pi{u) and Qa(u) are defined in equation ( |TT]) and 
equation (|16|) . 

(a) Since 

2(r 2 -n)(l +r 1 r 2 u 2 ) j u 2 k 2 \ j 2u k 2 \ 

[1 - (r 2 - n)u - r 1 r 2 u z \ l { k 3 h ) { k 3 k 3 J 

we have 

R' A (0) = 2(r 2 - n) - ^ = 2(r 2 - n) - - = 2(r 2 - n) - = 2(r 2 - n) f 1 - ^) . 

(17) 

Because the total mass of the belt is assumed to be less than the mass of the central star, 
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i.e. < M h < 1, we have R' A (0) < 0. Since R A (0) = 0, R' A (0) < 0, and R A ((l/r 2 ) _ ) > 0, 
we have that there exists -u* G (0, l/r 2 ) such that R A (u*) = 0, that is, Pj(w*) = Qa(w*).D 

(b) Since P^ ((l/?"2) + ) > and Pa((1/Vi)~) < 0; we have that there exists 
«*♦ G (l/r 2 , 1/ri) such that R A (u^) = 0, that is, = Qa^**)- 111 

(c) In the region of w > 1/ri, we have R A (l/ri) = Pj(\/r\) — exp j— ^ + f|^i} < 0. 

Further, because exp j— |^ + f|wj — ► as u — > oo, P/(w) — * 1 as u — > oo and R A (u) — > 1 
as m — > oo, there is a w*** > 1/n such that R A (u^) = 0, i.e. P/(it***) = <5a(«***)- d 

Figure l(a)-(b) are the numerical results for the solutions of equation flT5p. That is, 
we find all possible fixed points u for any given /c 2 and Mb with fixed values of r± and r 2 . 
Figure 1(a) are the results on k 2 — u plane, where dashed lines are for Mb = 0.1, dotted 
lines are for Mb = 0.5 and solid lines are for Mb = 1.0. There are three lines for each value 
of Mf, but only two solid lines can be seen because one solid line overlaps with the A; 2 -axis. 
Thus, there are three fixed points for any given values of A; 2 and M&. Figure 1(b) are the 
results on Mb — u plane, where dashed lines are for fc 2 = 0.1, dotted lines are for fc 2 = 0.5 
and solid lines are for fc 2 = 1.0. 

From these results, we can completely determine the locations of fixed points for 
different values of fc 2 and M&, where we assume < A; 2 < 1 and < M& < 1. These results 
are consistent with the analytic results of Theorem 3.1. 

3.2. Phase-Planes for Model A 

Following the linearization analysis, the eigenvalues A corresponding to the fixed points 
(u, v) satisfy the following equation 
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A 2 -(-l-^)=0, 



where A 1 = (In P/ («))/(«), so 



du 



lnPj(u) 

■u 2 



+ 



_L <9P/(m) 

uPiiu) du 



(18) 



Hence 



A = ±W-l-fc 3 



dA 1 



(19) 



9m 



We thus have two cases in the following: 

Case 1: If —1 — k 3 ^- > 0, then A are real with opposite sign, so in this case, the fixed 
point is a unstable saddle point. 

Case 2 : If —1 — k^^- < 0, then A are pure complex numbers with opposite sign, so in 
this case, the fixed point is a center point. 

Although we understand that when the real parts of the eigenvalues of a fixed point 
equal to zero in the first-order linearization analysis, the properties of this fixed point should 
be determined by the higher order analysis in general. However, to simply our language, we 
still use the term, center point, for any fixed point with zero real parts of eigenvalues. This 
is a good choice because our numerical results show that these points are in fact center 
points. 

In the following theorem, we discuss the properties of the fixed point for the physical 
region, i.e. u > l/r\. 

Theorem 3.2 

If the fixed point u* > 1/ri, then this fixed point u* is a center point. 



Proof : 
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First, from equation flI5D and equation flT8|), we calculate 



— 1 — k. 



dAi 
du 



-1 + k 



lnPj(u) 2k 3 (r 2 — Ti)(l + r\r 2 u 2 



-l + h 



u 2 u(r 2 u 2 — l)(r 2 u 2 — 1) 

k 2 l\ 2k 3 (r 2 - ri)(l + r x r 2 u 2 ) 



h + k 3 u 



u(r 2 u 2 — l)(r 2 u 2 



-2 + — 2k ^ T2 - r i)( 1 + r i r 2« 2 ) 



u(r 2 u 2 — l)(r 2 u 2 — 1) 



(20) 



For the convenience, we separate the proof into two parts: 

(i) when k 2 /2 < l/r\\ 
Since u* > l/r x , we have u* > k 2 /2. Thus, -2 + k 2 /u* < 0, so -1 - k 3 ^ 
fixed point u* is therefore a center point. 



< 0. The 



u=u* 



(ii) when k 2 /2 > l/r\\ 

Since u* > l/r 1 is a fixed point, u* satisfy Pi(u*) = Qa{u*). Since we consider the region 
u > 1/ri, we have 

< Pj(w) < 1 for u > 1/n (21) 

Moreover, because Qa only has one critical point at u — k 2 /2, i.e. Q' A {k 2 /2) = and 
Q" A {k 2 /2) < 0, QA{k 2 /2) is a global maximum and 

Qa (y) > Q A (k 2 ) = 1. 

Since Qa(0) = 1, QA{k 2 /2) > 1 and Q" A {u) < for u < k 2 , we have Qa(u) > 1 for u < k 2 . 
From equation (PH), we have RA(k 2 ) = Pi{k 2 ) — (^,4(^2) < and lim^oo Ra(u) = 1. 
Therefore, we have > k 2 . 

Because u* > k 2 /2, we have —2 + & 2 /^* < and thus —1 — k 3 ^- < 0. Therefore, the 
fixed point u* is a center point. □ 

The solution curves on the u — v phase plane are shown in Figure 1(c)- (d), where we 
set k 2 = 0.2, Mb = 0.3. In Figure 1(c), there are two vertical dotted lines, the left one is 
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u = l/r 2 and the right one is u — 1/ri, which divide the u — v plane into three regions. It 
is obvious that there is one fixed point in each region, which is consistent with Theorem 
3.1 and the fixed point in the region of u > \ jr\ is a center point, which is precisely what 
we have proved in Theorem 3.2. Figure 1(d) is just the detail of the solution curves near 
this fixed point in the region of u > l/r\. These figures reconfirm the analytic results and 
also make the behavior of the solution curve clear. We thus have known the behavior of the 
solutions completely. 



4. Model B 

In this model, we assume that the star-planet distance is larger than the outer radius 
of the belt, i.e. r > r 2 . Thus, / = f s + where f s and are defined in Section 2. ( f 8 
is defined by equation (H) and is defined by equation (|]). ) 

In this model, since r > r 2 > r 1; from equation (Eh in Section 2, we have 



fb 



TxGcm 



In 



(fx — r)(r + r 2 ) 



(r 2 — r)(r + r x ) 

By u—l/r and equation ([!]), we have the equation of motion: 

(1 - ri«)(l + r 2 u) 
(1 - r 2 u)(l + rxu) 



d 2 u 



-u 



Gm 2 TxGcm 2 f 
+ — Un 



Pu 



(22) 



Similarly, we set k 2 = Gm 2 /(l 2 ) and = (cirGm 2 ) / '(I 2 ) and transform equation 
to the following problem: 

du 

le = v 

^ = -u + k2 + -HPi(y)), 

do u 



(23) 
(24) 



where Pi(u) is defined in equation (|TT 
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4.1. Fixed Points for Model B 

The fixed points (u, v) of problem (Q) satisfy the following equations 



v = 



and 



u + k 2 + — ln(P/(u)) = 0. 



(25) 



u 



If we define 



Qb{u) = exp 




(26) 



and 



R b {u)=P i {u)-Qb{u), 



(27) 



then fixed points (u,Q) should satisfy Rb(u) = 0. In the following theorem, we discuss the 
properties of fixed points in different parameter space. 

Theorem 4.1 

For convenience, we define three regions as: 

Region (I): u < l/r 2 (this is the physical region for Model B), 

Region (II): l/r 2 < u < l/rj., 

Region (III): u > 1/ri. 

(a) In Region (I), if k 2 > l/r 2 , then there is no fixed point in this u < l/r 2 region. 

(b) In Region (II), there is at least one fixed point. 

(c) In Region (III), if k 2 < 1/Vi, then there is no fixed point. 
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(d) Within Region (III), if 1/n < k 2 /2 and Pi(k 2 /2) > Q B {k 2 j2), then there is one 
unique fixed point: u 2 E (1/n, k 2 /2) and at least one fixed point u 3 E {k 2 /2, k2). 

(e) Within Region (III), if 1/n < k 2 /2 and P/(A; 2 /2) < Q B (k 2 /2), then there is no fixed 
point in the region of (1/n , k 2 /2). 

(f) In Region (III), if 1/n < k 2 /2 and P/(u) > Q' B {u) for it e (A; 2 /2, k 2 ), then there is 
no fixed point. 

Proof : 

From the definition of Qb(u) in equation fl2Tf), we have Q' B {k 2 /2) = and Q" B {u) > 
for all u. Thus, Q B only has one critical point at k 2 /2 and Q B {k 2 /2) is a global minimum. 

Moreover, because Q B (0) = 1 = QB(k 2 ), we have 

Q B (u) < 1 for we (0,fc 2 ); (28) 
Q B {u) > 1 for u G (A;2 ? oo). (29) 

(a) Since Pr(0) = and P/(w) > 0, we have 

P z (it) > 1 for u E (o, — ^ . (30) 

Because of the condition fc 2 > l/ r 2) from equation (|2q) , we have 

Q B (u) < 1 for u E ^0, ^ . (31) 

From equation (|30D and (|31~D , we have R B [u) = Pi(u) — Q B {u) > for it < l/r 2 . 



Therefore, there is no fixed point in this physical region, u < l/r 2 .D 

(b) Because P j B((V r 2) + ) > and R B {l/ri) = — exp{rf/k 3 — k 2 /k 3 } < 0, there is at 
least one ui E (l/r 2 , 1/n) such that R B {u\) = 0. □ 
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(c) Because of the condition k 2 < 1/t 1; from equation (2S), we have 

Q B {u) > 1 for u G ,00) . (32) 

Since P/(l/ri) = 0, -P/(w) > and lim^oo Pi{u) = 1 , we have 

< Pr(u) < 1 for u G (— , 00) . (33) 



From equation fl32|) and (0), we have Rb(u) = Pi(u) —Qb( u ) < for all w G (1/ri, 00) 



Therefore, there is no fixed point in this case. □ 

(d) In the region of (l/r 1; k 2 /2), we have R B (l/n) < 0, R' B (u) = P'^u) - Q' B (u) > 
and from our assumption R B (k2/2) = Pj(k 2 /2) — QB(k 2 /2) > 0, so there is one unique fixed 
point u 2 G (l/ri,fc 2 /2). Because RB{k 2 ) < and Rs(k 2 /2) > 0, there is at least one fixed 
point u 3 G (k 2 /2, k 2 ). □ 

(e) Because R B {l/n) < 0, R' B {u) = P'^u) - Q' B {u) > for u G (1/^,^/2) and 
from our assumption RB(k 2 /2) = Pi(k 2 /2) — QB(k 2 /2) < 0, there is no fixed point in this 
(1/ri, k 2 /2) region. □ 

(f) First, we consider the region of (1/Vi, ^2)- Because of the condition 1/ri < /C2/2 and 
R' B (u) = Pj{u) — Q' B (u) > for u G (k 2 /2, k 2 ), from our assumption, we have 

= p;( M ) - Q' B (u) > for all u G f — , k 2 ) . 



Further, since i?s(l/ri) < and Rs(k 2 ) < 0, there is no fixed point in (1/n, k 2 ). 
Secondly, we consider the region of (fc 2 , 00). Because Pi{u) < 1, Qb{u) > Qs(k 2 ) = 1 and 
thus Rb(u) = Pi(u) — Qb( u ) < for all u > k 2 , so there is no fixed point in (fc 2 , 00) region. 
Therefore, from the above two, we have proved that there is no fixed point in the region of 
(l/ri,oo). □ 

In Theorem 4.1(a), we did not discuss the case when k 2 < l/r 2 . In Figure 2(a), we 
plot Rb as function of u for the cases that Mb = 0.5 and k 2 = 0.02 (solid line), k 2 = 0.04 
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(dotted line), k 2 = 0.06 (dashed line), k 2 = 0.08 (long dashed line). All these values of k 2 
satisfy k 2 < l/r 2 . It is obvious that there is no root for k 2 = 0.08 (long dashed line) and 
two roots for others. From these numerical results, we know that there could be two fixed 
points or no fixed point in the (0, l/r 2 ) region if k 2 < l/r 2 . 

On the other hand, in Theorem 4.1(c)-(f), we did not discuss the case when 
k 2 /2 < 1/ri < k 2 . In Figure 2(b), we plot Rb as function of u for the cases that Mb = 0.5 
and k 2 = 0.4 (solid line), k 2 = 0.5 (dotted line), k 2 = 0.6 (dashed line). All these values of 
k 2 satisfy k 2 /2 < 1/ri < k 2 . It is obvious that there is no root for k 2 = 0.4 (solid line) and 
two roots for others. From these numerical results, we found that there could be two fixed 
points or no fixed point in in Region (III), i.e. (1/rb oo) if k 2 /2 < 1/n < k 2 . 

Figure 2(c)-(d) are the numerical results for the solutions of equation (p5|) . That is, we 
find all possible fixed points u for any given k 2 and M& with fixed values of r\ and r 2 . 

Figure 2(c) are the results on k 2 — u plane, where dashed lines are for Mb = 0.1, dotted 
lines are for Mf, = 0.5 and solid lines are for M& = 1.0. From this figure, we know that there 
are three fixed points when k 2 is closer to or 1, i.e. the left or right side of the figure. 
However, there could be one fixed point only in the middle. Figure 2(d) are the results on 
Mb — u plane, where dashed lines are for k 2 = 0.1, dotted lines are for k 2 = 0.5 and solid 
lines are for k 2 = 1.0. 

Checking the detail of these figures and Theorem 4.1, we found that they are completely 
consistent with each other. 

4.2. Phase-Planes for Model B 

To study the phase planes, we need to know the eigenvalues of every fixed point. The 
eigenvalues A corresponding to a fixed point (u, v) satisfy the following equation: 
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A 2 -(-l + ^) = 0, 

where ?* = *Mfi(«))/tt] = "jn^u) + 1 d { P l{ u)) 
ou ou u uPAu) ou 



Hence 



We have two cases for the eigenvalues in the following: 

Case 1: If — 1 + ^3^7 > 0, then A are real with opposite sign, so the fixed point is a 
unstable saddle point. 

Case 2: If —1 + ks^- < 0, then A are pure complex numbers with opposite sign, so the 
fixed point is a center point. 

Theorem 4.2 

If k 2 r 2 < 2, then there is a unique fixed point U\ G (l/r 2 , 1/ri) and this fixed point is a 
center point. 

Proof : 



Since fixed point (u,v) satisfy equation ([25]), from equation fl34|), we have 

_ 1 + h dB 1 = J InPjju) | 2 (r 2 -r 1 )(l + r 1 r 2M 2 ) | 

<9w \ u 2 u (1 — r 2 -u 2 )(l — r 2 u 2 ) J 

+ f 1 fc 2 2 (r 2 -ri)(l + rir 2 n 2 ) | 

3 \ ^3 k 3 u u (1 — r 2 w 2 )(l — r 2 u 2 ) J 
_ 2 ^2 2fc 3 (r 2 - ri)(l + r x r 2 u 2 ) 
u u (1 — r%u 2 )(l — rfu 2 ) 

^ 2 k 2 M b (l + r\r 2 u 2 ) , , 

= -2 + -- , 2 tttt 2 2V (36) 

Since fc 2 /2 is a global minimum of Qb{u) and A; 2 /2 < l/r 2 , we know that Q'b{u) > 
for u G (l/r 2 ,l/ri). Because P'j(u) < and Q' B {u) > for all u G (l/r 2 ,l/ri), we have 
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R' B {u) = P'j(u) — Q' B (u) < for all u G (l/r 2 ,l/r 1 ). Moreover, it can be shown that 
-Rb((1/t 2 ) + ) = oo and R B {l/ri) = —Q B {\/ri) < 0, thus there is a unique fixed point 
U\ G (l/ r 2j V r i) sucn ^at R B {u\) = 0. 

Because this fixed point U\ G (l/r 2 , 1/Vi), we have 

-2 + — < -2 + k 2 r 2 < 0. 

Ml 



Since the third term of equation (36) is always negative, we have —1 + k 3 ^- 



< 0. 

U=U\ 



Therefore, this fixed point is a center point. □ 



The solution curves on the u — v phase plane are shown in Figure 3(a)-(d), where the 
vertical dotted lines are u = l/r 2 and u = 1/t\. Figure 3(a) is the result when k 2 = 0.05, 
Mb = 0.5, which corresponds to the left region of Figure 2(c) (on the dotted line and k 2 
is closer to ). Thus, there are three fixed points: one is larger than l/r 2 and others are 
less than l/r 2 . Figure 3(b) is the result when k 2 = 0.8, M b = 0.5, which corresponds to the 
right region of Figure 2(c) (on the dotted line and k 2 is closer to 1 ). Thus, there are three 
fixed points: one is less than 1/n and others are larger than l/r\. Figure 3(c) is the result 
when k 2 = 0.5, Mb = 0.5, which corresponds to the middle region of Figure 2(c) (on the 
dotted line and about the region between 0.08 and 0.55) Thus, there is only one fixed point. 
Figure 3(d) is the result when k 2 = 1.0, = 0.8, which corresponds to the right region of 
solid lines in Figure 2(d). Since is closer to 1, there are three fixed points: one is less 
than 1/ri and others are larger than \jr\. Therefore, the numbers of fixed points in Figure 
3 are in fact completely consistent with Figure 2. 

Figure 3(a), 3(c) and 3(b) give us a chance to see the transition of phase plane along 
the dotted line in Figure 2(c) from left side to the right side. In Figure 3(a), there are 
one center and one saddle point in the region of u < l/r 2 and one center in the region of 
l/r 2 < u < \jr\. However, in Figure 3(c), the left two fixed points disappear, so there is 
only one fixed point, which is the center in the region of l/r 2 < u < 1/ri. Moreover, the 
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shape of the solution curves in the region of u > l/r% already become different from the 
curves close to a center. Finally, in Figure 3(b), there is still one center point in the region 
of I/V2 < u < 1/ri but there are two new fixed points: one saddle and one center in the 
region of u > l/r\. The locations of these bifurcation points, where transitions of phase 
plane occur, can be easily seen from Figure 2(c). 

These figures reconfirm the analytic results and also make the behavior of the solution 
curves clear. 



5. Model C 



In this model, we assume that the star-planet distance is between the inner and outer 
radius of the belt, i.e r\ < r < r 2 . Thus, / = f s + f a + fb e where / s , f a and /& e are defined 
in Section 2. 



From equation (|l|), set u = 1/r, we have 



de 2 



-u + 



Gm 2 ac du TcGcm 2 



I 2 lu d9 



l 2 u 



In 



eV 



(1 + nu)(l + r 2 u) 
(1 — r\u)(r2U — 1) — e 2 u 2 



(37) 



To simplify the equation, we define 

(1 + r\u)(l + r 2 u) ( e 2 u 2 



PAu) = 



[1 — riu){r2U — 1) \4 — e 2 u 2 



(38) 



Hence, equation ( p7[ ) becomes 

d 2 u k\ du 
dO 2 + ~^d8 



-u + k 2 -—\n(P e (u)) 
u 



(39) 



where k\ = (ac)//, k 2 = (Gm 2 )/(l 2 ) and = (c7iGm 2 )(l 2 ). Therefore, we transform 
equation fl3"9p to the following problem 
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du 

Te = 

= _ U + / C2 _^_^l n p e w . 40 
du u u 



5.1. Fixed Points for Model C 



The fixed points (u,v) of problem (40) satisfy the following equations 



We define 



and 



v = 

and -u + k 2 - — v- — In(PJu)) = 0. (41) 
it u 



g c (tt) = exp|-^ + ^«|, (42) 



it: c («) = p e ( w )-g c («), (43) 



then the fixed points (w, u) satisfy R c (u) = 0. In the following theorem, we discuss some 
properties about the fixed points. 

Theorem 5.1 

(a) There is at least one fixed point u\ G (0, 1/7*2) ■ 

(b) There is at least one fixed point u e * G (l/r 2 ,u), where u G (1/7*2, V r i) suc h that 
P e '(7i) = 0. Moreover, 

m £ * = h 0(e 2 ), so limu e * = — . (44) 

r 2 e-^o r 2 



-21 - 



(c) There is at least one fixed point u* G (u, 1/n), where u G (l/r 2 , 1/ri) such that 
P' e {u) = 0. Moreover, 

< = — + 0(e 2 ), so limit* = — . (45) 

(d) There are at least two fixed point u 2 ,u 3 G (l/ri,2/e). 
Proof : 

(a) Since P c (0) = — Q c (0) < and P c ((l/r 2 ) _ ) > 0, there is a root U\ G (0, l/r 2 ) such 
that R c {ui) = 0. Thus, tti is a fixed point. □ 

(b) From equation fl38|), if we consider w G (l/r 2 , l/?"i), then we have 

p// n _ 2e 2 M f (n + r 2 )(l - r 1 r 2 u 2 )u 4(1 + nw)(l + r 2 u) 1 



(1 — r 1 u)(r 2 u — 1)(4 — e% 2 ) | (1 — r 1 -u)(r 2 M — 1) (4 — e 2 w 2 ) 



(46) 



Since P e '((l/r 2 )+) < 0, P e '((l/n)-) > and P e "(«) > for u G (l/r 2 , 1/n), there is a 
w G (l/r 2 , 1/n) such that P £ '(w) = 0. 



Since P £ (tt) = 0, from equation (|46|) , we have 

(fi + r 2 )(l - r 1 r 2 (u) 2 )u _ 4(1 + n«)(l + r 2u) 
{1 - r x u){r 2 u - 1) (4-e 2 (w) 2 ) 

Therefore, from equation (f47|) , we have 

(l + nM)(l + r 2 M) / 6 2 (^) 2 \ 
el ; {l-r x u){r 2 u- 1) ^4-e 2 (w) 2 / 
4e 2 w(l + ri«) 2 (l + r 2 w) 2 



(47) 



(4£ 



(4 — e 2 (-u) 2 ) 2 (r! + r 2 )(l — rir 2 (u) 2 ) 

From equation fl4"8|), we have P e (u) = 0(e 2 ). 
Since P c ((l/r 2 )+) > 0, and 

R c (u) = P e {u) - Q c (u) = 0(e 2 ) - Q c (u) < 0, (49) 
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we have that there is at least one fixed point u e * G (l/r 2 ,u) such that R c (u t *) = 0. 



In the following, we prove that equation ( fPj) holds. Since R c {u e *) = 0, we have 
P e (u € *) = Q c (u e *), that is, 

J(l + riu e *)(l+r 2 u e *) / e 2 u% \\ 1 2 k 2 
'"{ (I - r lM „)(r 2t .„ - 1) IrAl j | = + (50) 

We found the right hand side in equation ([5(]) is bounded for all e > 0, since 
■u e * G (l/r 2 ,-u). We can rewrite the left hand side of equation (|50D to be 

1 1 - riu e * J \ 4 - e 2 u 2 * J \ r 2 u e * - 1 J 
It is easy to show that the first term and the second term in equation fl5"T|) are bounded for 
all e > since u e * G (l/r 2 ,tt). Therefore, the third term in equation flST]) is bounded. That 
is, there exists a M > such that 

- M < ln { e2M " 1 < M. (52) 
I r 2 w e * - 1 I 



Since u e * G (l/r 2 ,w), from equation (p2|), we have 



1 + /-, /-\9 6 / Ha < r 2 M « < 1 + 6 JTTs • (53) 



;i/«) 2 exp(-M) exp(M) 



Therefore, w e * = ^ + 0(e 2 ) and 

limw e * = — .□ 

e^O r 2 



c) From equation ([49]), we have R c (u) < and i? c ((l/ri) ) > 0, so there is at least 



one fixed point u* G (u, 1/ri) such that R c (u*) = 0. 

For equation (f|^), we have the same argument as the proof of (b).D 

(d) From equation fl38|), if we consider w G (1/ti, 2/e) then we have 

2e 2 w J (r 1 + r 2 )(r 1 r 2 u 2 )u - 1 4(1 + tyu)(1 + r 2 w) 1 

e ^ ~ (r lM -l)(r 2 M-l)(4-e% 2 ) \ (^-1)^-1) (4 - e 2 w 2 ) J ' 



(54) 
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Since P^l/r^) < 0, P e '({2/e} _ ) > and P e "(u) > for w G (l/r 2 , 1/n), there is a 
u G (l/r 1; 2/e) such that P' e (u) = and P e (w) is a minimum value for u G (l/ri,2/e). 



Since P e '(w) = 0, from equation (|38|), we have 

(ri + r 2 )(rir 2 (-u) 2 - l)u _ 4(1 + riw)(l + r 2 w) 
(nri- l)(r 2 w- 1) ~ (4 - e 2 (u) 2 ) 

Therefore, from equation fl55|), we have 



(55) 



= (l + nfl)(l + r 2 fl) / e 2 (u) 2 \ 
eK ' " (nu - l){r 2 u - 1) V 4 - e2 W 2 / 
4e 2 u(l + 7Yu) 2 (l + r 2 u) 2 
(4 — e 2 (u) 2 ) 2 (ri + r 2 )(rir 2 (u) 2 — 1) ' 



(56) 



and we found the minimum value P e (u) = 0(e 2 ). Therefore, R c (u) = P e {u) ~ Q(u) < 0. 
Since -R c ((l/r 1 ) + ) > 0, P c ((2/e)~) > and R c {u) < 0, there are at least two root 
u 2 G (l/r^u) and u 3 G (u,2/e) such that R c (u 2 ) = R c (u 3 ) = 0. □ 



5.2. Phase-Planes for Model C 



Next, we find the eigenvalues A satisfy the following equation 

dC x . 



ac , 



A + —A - (-1 - k 3 — 

IU OU 



0. 



(57) 



Hence, 



where C\ 
dd _ d[\n(P e (u))/u] 



ln(P e (w)) 



du du 
Thus, from equation (|57|) , we have 



u 



lnPe(«) 1 d{P e {u)) 



uPJu) du 



fau\ 2 ( , dC x 



du 



(51 



(59) 



(60) 
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Let A = 1 — 4 (l + ^3^7) (cui) ' then we have three possible cases : 

Case 1: If A < 0, then two eigenvalues Aj.,2 = —a ± bi for a, b > 0, so the fixed point is 
a stable spiral point. 

Case 2: If < A < 1, then two eigenvalues are negative, so the fixed point is a stable 
node. 

Case 3: If A > 1, then one eigenvalue is positive and the other is negative, so the fixed 
point is a saddle point. 

Theorem 5.2 

(a) The fixed point, u e * = ^- + 0(e 2 ) is a saddle point. 

(b) The fixed point, u* = j- + 0(e 2 ) is a stable spiral point. 
Proof : 

Since fixed point (u,v) satisfy equation (PD , from equation (|59D, we have 
dd 



du 



J lnP e (w) 2(ri + r 2 )(l - r 1 r 2 u 2 ) 

1 + ft 3 \ n 7^ o nw o n 7T + 



w 2 u(l — r 2 u 2 )(r 2 u 2 — 1) u 2 (A — e 2 u 2 )j 
_ 2 ^2 2fc 3 (ri + r 2 )(l -rir 2 w 2 ) | 8A; 3 

K tt(l — r 2 u 2 )(r 2 u 2 — 1) u 2 (A — e 2 u 2 ) 

In order to know which case it should be for different fixed point, we further calculate: 

\ ou ) \ac J 

. ( 9 _ ^2 _ 2fc 3 (n + r 2 )(l - rir 2 w 2 ) 8A; 3 \ ( lu\ 2 



11 

2 



■u(l — r 2 u 2 )(r 2 u 2 — 1) m 2 (4 — e 2 w 2 ) y \oic / 



lu\ f4k 2 l 2 u\ 8k 3 l 2 u( y r 1 +r 2 )(l-r 1 r 2 u 2 ) 32k 3 l 



2 



■ (62) 



\ac/ \ a 2 c 2 / a 2 c 2 {l — r 2 u 2 ){r 2 u 2 — 1) a 2 c 2 (4 — e 2 w 
(a) We consider the fixed point, tt e * = + 0(e 2 ). Since the fourth term of equation 
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(|62|) is positive with magnitude 0(l/e 2 ) and other terms are bounded , thus A > 1 and the 
fixed point u e * = i + 0(e 2 ) is a saddle point. □ 

(b) We consider the fixed point, u* = y- + 0(e 2 ). Since the fourth term of equation 
(|62|) is negative with magnitude 0(l/e 2 ) and other terms are bounded , thus A < and the 
fixed point u* = + 0(e 2 ) is a stable spiral point. □ 

The solution curves on the u — v phase plane are shown in Figure 4(a)-(d), where we 
set ki = 0.1, M b = 0.5 and e = 10~ 3 . In Figure 4(a), there are two vertical dotted lines, the 
left one is u = 1 jr-i and the right one is u — 1/ri , which divide the u — v plane into three 
regions. Figure 4(b) is the detail for the region near the line of u = I/V2 in Figure 4(a), 
which make it possible to see the fine structures. We can see that there is one center and 
one saddle point near the dotted line u — l/r 2 . Figure 4(c) is the detail for the region near 
the line of u — 1/ri in Figure 4(a). We can see that there is one spiral point and one saddle 
point near the dotted line u = l/r\. Figure 4(d) is a closer-looking of the spiral point in the 
region of u > \jr\. 

These figures reconfirm the analytic results and also make the behavior of the solution 
curve clear. 

6. Conclusions 

We have studied the properties of solution curves on the phase plane of dynamical 
systems of planet-belt interaction by the standard fixed-point analysis. 

The system is divided into three regions: (a) the region between the central star and 
the inner edge of the belt, (b) the region out of the outer edge of the belt, (c) the region 
between the inner edge and the outer edge of the belt. The dynamics in three regions are 
governed by three different models and the planet moves around these regions. 
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We analytically prove some properties of the system and also show some bifurcation 
diagrams of fixed points for different cases: the number and properties of fixed points 
changed with the values of parameters k 2 and M& in general. 

Our analytical and numerical results show that, in most cases, the locations of fixed 
points depend on the values of these parameters and these fixed points are either center 
(structurally stable) or saddle (unstable). 

On the other hand, we found that there are two special fixed points: one is on the 
inner and another is on the outer edges of the belt. The one on the inner edge is a stable 
focus (asymptotically stable) and the one on the outer edge is a saddle (unstable). 

This interesting result is consistent with the observational picture of Asteroid Belt 
between the Mars and Jupiter: the Mars is moving stablely close to the inner edge but the 
Jupiter is quite far from the outer edge. 
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Fig. 1. — Bifurcation diagrams of fixed points and phase planes for Model A with n = 3, 
r 2 = 6. (a) The bifurcation diagram of fixed points on k 2 — u plane, where dashed lines 
are for M b = 0.1, dotted lines are for M b = 0.5 and solid lines are for M b = 1.0. (b) The 
bifurcation diagram of fixed points on M b — u plane, where dashed lines are for k 2 = 0.1, 
dotted lines are for k 2 = 0.5 and solid lines are for k 2 = 1.0. (c) The solution curves on the 
u — v plane, where we set k 2 = 0.2, M b = 0.3. There are two vertical dotted lines, the left 
one is u = l/r 2 and the right one is u — \ jr\. (d) The detail near the rightest fixed points 
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Fig. 2. — Bifurcation diagrams of fixed points for Model B with r\ = 3, r 2 = 6. (a) R B 
as a function of u for the cases that M b = 0.5 and k 2 = 0.02 (solid line), k 2 = 0.04 (dotted 
line), k 2 = 0.06 (dashed line), k 2 = 0.08 (long dashed line), (b) R B as a function of u for the 
cases that M b = 0.5 and k 2 = 0.4 (solid line), k 2 = 0.5 (dotted line), k 2 = 0.6 (dashed line), 
(c) The fixed points on k 2 — u plane, where dashed lines are for M b = 0.1, dotted lines are 
for M b = 0.5 and solid lines are for M b = 1.0. (d) The fixed points on M b — u plane, where 
dashed lines are for k 2 — 0.1, dotted lines are for k 2 = 0.5 and solid lines are for k 2 = 1.0. 
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Fig. 3. — The phase planes for Model B with r\ = 3, r 2 = 6. (a) The solution curves on the 
u — v plane when k 2 — 0.05, M b = 0.5. The vertical dotted line is u = l/r 2 . (b) The solution 
curves on the u — v plane when k 2 = 0.8, M b = 0.5. The left vertical dotted line is u — l/r 2 
and the right vertical dotted line is u — l/r±. (c) The solution curves on the u — v plane 
when k 2 = 0.5, M b = 0.5. The left vertical dotted line is u — l/r 2 and the right vertical 
dotted line is u — l/r\. (d) The solution curves on the u — v plane when k 2 = 1.0, M b = 0.8. 
The left vertical dotted line is u = l/r 2 and the right vertical dotted line is u — \ jr\. 
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Fig. 4. — The phase planes for Model C when T\ = 3, r 2 = 6, k 2 = 0.1, M b = 0.5 and 
e = 10~ 3 . (a) The u — v phase plane, the left vertical line is u — l/r 2 and the right one is 
u — \jr\. (b) The detail for the region near the line of u — l/r 2 in (a), (c) The detail for 
the region near the line of u — l/r± in (a), (d) A closer-looking of the focus point in the 
region of u > l/r\. 



